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Abstract 

We make a detailed comparison of two tadpole renormalization schemes in 
the context of the quarkonium hyperfine splittings in lattice NRQCD. We 
renormalize improved gauge-field and NRQCD actions using the mean-link 
«o,l in Landau gauge, and using the fourth root of the average plaquette 
ito,p- Simulations are done for the three quarkonium systems cc, be, and bb. 
The hyperfine splittings are computed both at leading (0(Mqv 4 )) and at 
next-to-leading (0(Mqw 6 )) order in the relativistic expansion, where Mq is 
the renormalized quark mass, and v 2 is the mean-squared velocity. Results 
are obtained at a large number of lattice spacings, in the range of about 
0.14 fm to 0.38 fm. A number of features emerge, all of which favor tadpole 
renormalization using uq l- This includes much better scaling behavior of 
the hyperfine splittings in the three quarkonium systems when uq^l is used. 
We also find that relativistic corrections to the spin splittings are smaller 
when uo t L is used, particularly for the cc and be systems. We also see signs 
of a breakdown in the NRQCD expansion when the bare quark mass falls 
below about one in lattice units. Simulations with uq^l also appear to be 
better behaved in this context: the bare quark masses turn out to be larger 
when tio,L is used, compared to when uo,p is used on lattices with comparable 
spacings. These results also demonstrate the need to go beyond tree-level 
tadpole improvement for precision simulations. 
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I. INTRODUCTION 



Tadpole improvement of lattice actions |TJ has become an essential ingredient in numerical 
simulations of hadronic systems. Tadpole improvement has revitalized interest in simulations 
on coarse lattices, and is playing an important role in current efforts to extract continuum 
results from simulations on fine lattices. 

Tadpole diagrams in lattice theories are induced by the nonlinear connection between 
the lattice link variables and the continuum gauge fields. This has been shown to cause 
large radiative corrections to many quantities in lattice theories. Fortunately, most of the 
effects of tadpoles can be removed by a simple mean field renormalization of the links |TJ 

um - (i) 

MO 

where an operator dominated by short- distance fluctuations is used to determine uq. 

One of the earliest applications of tadpole improvement was in the development of lattice 
nonrelativistic quantum chromodynamics (NRQCD) 0-0]. Precision simulations of the T 
system in NRQCD have provided many important phenomenological results, including the 
strong coupling constant ||] and the 6-quark pole mass |5]||. However the situation for 
charmonium is more problematic, due to large relativistic corrections in this system Q. 

In fact the quarkonium spectrum provides a powerful probe of tadpole renormalization. 
The quarkonium fine and hyperfine spin splittings in particular are very sensitive to the 
details of the NRQCD Hamiltonian, with the relevant operators undergoing large tadpole 
renormalizations. For example, it has been shown || that scaling of the charmonium hy- 
perfine splitting is significantly improved when the tadpole renormalization is determined 
using the mean-link uq l measured in Landau gauge M : 



W0)i = (|ReTrC/ M ), 0^ = 0, (2) 
compared to when the fourth root of the average plaquette m ,p is used: 

u ,P = (iReTrf/ pl ) 1/4 . (3) 

Equation (^) has been employed in most previous lattice simulations. First evidence for a 
more continuum- like behavior of lattice actions using m 0j l came from studies of rotational 
symmetry restoration in the heavy quark potential 0. More recently, improved scaling of 
the charmonium spectrum from relativistic actions HlPl , and of the SU(2) glueball spectrum 



nj, have also been observed in simulations using uo.z,. 

In this paper we make a detailed comparison of the two tadpole renormalization schemes 
mo,l and uo,p (when implemented at tree-level) in the context of the quarkonium hyper- 
fine splittings in NRQCD. This is done for the three quarkonium systems cc, be, and bb. 
The hyperfine splittings are computed both at leading {0(Mqv a )) and at next-to- leading 
(0(Mqv 6 )) order in the relativistic expansion, where Mq is the renormalized quark mass, 
and v 2 is the mean-squared velocity. Results are obtained at a large number of lattice spac- 
ings, in the range of about 0.14 fm to 0.38 fm. All quantities are calculated in the two 
tadpole schemes after careful re-tuning of the lattice action parameters for each system. 
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We find that scaling of the hyperfine splittings in all three quarkonium systems is signif- 
icantly improved when m 0i l is used. This confirms and extends the scaling analysis done in 
Ref . |J . In addition we see signs of a possible breakdown of the NRQCD effective action at 
smaller lattice spacings, when the bare quark mass in lattice units aMg falls below about 
one (scaling studies of the T spectrum, as a probe of the NRQCD expansion, have previ- 



ously been reported in Ref. |12fl). Large changes in the cc and be splittings with mq,p are 



evident at smaller lattices spacings. Once again, simulations with m ,l are better behaved 
in this context: the bare charm-quark mass turns out to be much larger when u 0j l is used, 
compared to when u 0j p is used on lattices with comparable spacings. Although there are 
clear pathologies in the iio,p data at small bare masses, our conclusions here are necessarily 
tentative, since we are unable to do realistic simulations on sufficiently fine lattices to get 
aM c ° appreciably below one with uo,l tadpole improvement. 

We also find that the choice of tadpole renormalization scheme has a strong effect on the 
apparent convergence of the velocity expansion which underlies the NRQCD effective action. 
We find that the relativistic corrections to the hyperfine splittings are smaller when m 0j l is 
used, particularly for the cc and be systems. This casts new light on the results obtained 
in Ref. ||, where relativistic corrections to spin splittings in NRQCD were first calculated. 
Results obtained here show that the charmonium hyperfine splitting is reduced by about 
30% when u 0i l is used to renormalize the action. This is consistent with a naive estimate of 
v 2 E from operator expectation values in NRQCD ||. Hence while the velocity expansion for 
charmonium is subject to large corrections, it may not be as unreliable as was suggested in 
Ref. P| , where simulations with uq } p on relatively fine lattices showed relativistic corrections 
of about 60% in the charmonium hyperfine splitting. (Relativistic corrections have more 
recently been analyzed in the T system ]13|,14|, and have also been studied in heavy-light 
mesons [|H)|-|T7|] ) . 

Although these results clearly favor using uq.l for tadpole improvement, they also serve 
to underline the need to go beyond tree-level matching of lattice actions. Results from 
this study and others (see e.g. Refs. IPHT^I) demonstrate that precision results can only be 
obtained once uncertainties due to 0(a s ) renormalizations are removed (for some recent 
work in this connection see e.g. Refs. [|I8|,|I9"|1 ) . 



II. DETAILS OF THE SIMULATIONS 



A. Lattice Actions 



The lattice NRQCD effective action for quarkonium is organized according to an ex- 
pansion in the mean squared velocity v 2 of the heavy quarks, with corrections included for 
lattice artifacts. The effective action, including spin-independent operators to 0(t> 4 ), and 
spin-dependent interactions to 0(t> 6 ), was derived in Ref. M. Following Refs. P,f7|, we use 
the evolution equation 

Gm=(l-^)V|(l-^) n (l-a^)G 4 (t>l), (4) 
where the initial evolution is set by 
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G^U-^YuUl-^) &,„. (5) 



2n J * V In 
On the lattice the leading kinetic energy operator Hq is given by 

A( 2 ) 

where Mq is the bare quark mass and A^ 2 ) is the lattice Laplacian. 

Relativistic corrections are organized in powers of the heavy quark velocity: 

5H = 5H (4) + 5H (6) . (7) 

contains spin-independent relativistic corrections and leading-order spin interactions: 

-ca ,? rn „ (T • (A x E — E x A) — c 4 -^<t ■ B 

8(M°) 2 1 } 4 2M° 

a 2 A( 4 ) a(A( 2 )) 2 
+ ° 5 24M0 C6 16n(M°) 2 ' ( } 

with the last two terms coming from finite lattice spacing corrections to the lattice Laplacian 
and the lattice time derivative respectively. The parameter n is introduced to remove insta- 
bilities in the heavy quark propagator caused by high momentum modes || . Spin-dependent 
relativistic corrections for quarkonium first appear at 0(v 6 ): 



-eg 



3g 



64(M°) 



{A (2) ,<r- (A x E-E x A)}, 



Spin-independent corrections at 0(v 6 ) are not considered here (these operators may in fact 
have indirect effects on spin splittings 0,0). As in Ref. || simulations were done with the 
derivative operators and the clover fields corrected for their leading discretization errors. 
This is indicated by the tilda superscripts on these operators in Eqs. (^) and (^j). Complete 
expressions for the operators can be found in Refs. |||8|]. 

At tree-level all of the coefficients q in Eqs. (§) and (|9]) are one. However tadpole 
improvement is crucial in order to eliminate large radiative corrections. Independent sets of 
simulations were done using the two tadpole renormalization schemes described in Section I: 
the mean-link in Landau gauge, Eq. (§) (using a standard lattice implementation of the 
continuum Landau gauge fixing J2UJ), and the fourth root of the average plaquette, Eq. (^[). 



The links were rescaled in the simulation before they were input to the quark propagator 
subroutine, to be sure that Eq. (|l|) was correctly implemented in all terms in the heavy quark 
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action. The gauge-field configurations were generated using an 0(a 4 )-accurate tadpole- 
improved action [21] 

S[U] = /?£ I^Tr (l - U pl ) - E l ReTr (! " ^rt) > (10) 
pi rt 

where the sums are over all oriented lxl plaquettes (pi) and 1x2 rectangles (rt). 

Meson creation operators were constructed from quark (ip^) and antiquark (x*) creation 
operators @,|||: 

^(x)T(x)x\x), (11) 



using a gauge-invariant smearing function ||22[ 

T(f) = 7 t (x)n(f) 7 (x), (12) 
where the 2x2 spin matrix Q(x) selects the quantum numbers of interest, and 

1 (x) = (l + eA^(x)Y a . (13) 



The weight e and the number of smearing iterations n s were adjusted to optimize the overlap 
with the ground state. 

Meson correlation functions ||[7j were computed for the 1 S'o (12 = I), 3 Si (Q = cii) 
and l P\ (Q = Aj) mesons. The three triplet P-wave correlators ( 3 Po, 3 Pi, 3 P2) were also 
analyzed, but this data is of insufficient quality to report here. Correlation functions were 
evaluated at both zero momentum and at the smallest allowed nonzero momentum. 



B. Simulation Parameters 

Six lattices were generated using the mean link in Landau gauge to set the tadpole factor 
(mo,l) and seven lattices with comparable spacings were generated using average plaquette 
tadpoles (ito,p)- The parameters of these thirteen lattices are given in Tables [I] and [TT]. Note 
that we use /3l to denote the lattice coupling for simulations with M i , and j3p for lattices 
with Mo.p- To check that finite volume effects are not an issue on the lattices with the 
smallest spacings, we performed some runs at (3p = 7.3 on a 16 4 volume, and found no 
significant change in the results. 

A standard Cabbibo-Marinari pseudo heat bath was used to generate the gauge field 
configurations. The number of updates between measurements varied from 10 for lattices 
with the largest spacings to 20 for the smallest; autocorrelation times satisfied r < 0.5 in 
all cases. Smeared-smeared correlators were used, with typically 5-10 smearing iterations, 
and a smearing weight e = 1/12 was used in all cases. 

The lattice spacings were determined from the spin-averaged IP — IS mass difference, 
which we set to 458 MeV (the experimental value for charmonium); this mass difference is 
thought to be about the same for all quarkonium systems ||23|| . The difference between the 
singlet 1 Pi and the spin-averaged 3 Si, x So states was used for this purpose. 
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The bare quark masses were tuned by calculating the kinetic masses M^ n of the 1 So 
states in physical units, which were extracted from fits to the energies Ep of the boosted 
states using the form 

E *- E ° = m^- (14) 

Fits were made to the state with momentum components (1,0,0) in units of 2n/(Na); 
relativistic corrections to the dispersion relation made little difference in the fit values of Mk; n . 
The bare quark masses Mq were fit to the following 1 Sq kinetic masses: M c5 = 2.98 GeV, 
M b5 = 6.28 GeV, and M h i = 9.46 GeV. The value of the be S'-wave meson mass adopted 
here was obtained in a previous NRQCD analysis of the be system ||23|| . 

A potential complication arises in that the lattice spacing and quark masses determined 
from different systems on the same quenched configurations do not generally agree. Since we 
compare results for a given system obtained over a wide range of lattice spacings, we have 
attempted to minimize systematic effects from tuning errors by re-tuning the parameters 
for each system on each ensemble of configurations. In the case of the be system we used 
the c-quark mass determined from the M c5 kinetic mass, and tuned the 6-quark mass to 
reproduce Mb 5 . Note that the 6-quark mass was re-tuned in simulations of the bb system, 
to obtain the correct value of M b i. The final simulation results for the kinetic masses are 
accurate to within 3% in all cases. 

The resulting quark masses and lattice spacings for the three quarkonium systems for 
the NRQCD action at 0(v 6 ) are given in Tables [III] and |IV]. Note again that the lattice 



spacing is given separately for each system, the differences possibly reflecting effects due to 
quenching, which has been conjectured to play a role in setting the scale in these and other 
hadronic systems (see e.g. Ref. ||23|| ). The 6-quark mass is similarly given separately for the 
be and bb systems. We found only small changes in the quark masses and lattice spacings 
when they are determined from the NRQCD action at 0(v 4 ). 

Although we think that it is worthwhile to minimize systematic errors by re-tuning the 
quenched lattice parameters for the different quarkonium systems, and have done so in all 
simulations reported here, these effects are actually small. It is important to note that none 
of the conclusions reached in this study are changed if our re-tuning procedure is modified. 



III. RESULTS AND ANALYSIS 

Hyperfine splittings for the three quarkonium systems cc, be, and bb were extracted on the 
six lattices with uq,l (cf. Table |), and the seven lattices with uo t p (cf. Table |I|). To illustrate 
the quality of the raw data we show effective mass plots m e g(T) = — \og{G(T) / G{T — 1)) for 
the be system at 0(v 6 ), on the two lattices with the smallest spacings, in Fig. [I] (j3 L = 7.5) 
and Fig. |(/? P = 7.3). 

Single exponential fits to the correlation functions were used to get the masses of the 
ground states, and a jackknife analysis was used to estimate statistical errors. The hyperfine 
splittings were obtained from a fit to the ratio of 3 Si and 1 Sq correlation functions. Detailed 
fit results corresponding to the data in Figs. |T] and [2] are given in Tables [V] and |VT|. Final 
estimates of the dimensionless energies were obtained by finding two or three successive 
t min /t max intervals for which the fit results overlap within statistical errors. The largest 
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statistical error in the overlapping fits was used as an estimate of the error in the final fit 
value. We give the final fit results for the dimensionless energies for the be states at 0(v 6 ) 
in Tables |VTJ and |VTTI| . 

The final fit results for all hyperfine splittings in physical units are given in Tables [IX 
and £X|, where the quoted errors are purely statistical. We expect a systematic error of order 
10% in the hyperfine splittings, coming from uncertainties in the quark mass determinations 
due to 0(v 6 ) spin-independent relativistic corrections • Effects due to quenching may 



also be large in the case of the charmonium hyperfine splitting [24], as the estimates from 



relativistic p^|TU| and NRQCD || actions are all substantially smaller than the experimental 
value of (118 ± 2) MeV. 

The splittings for the cc, be, and the bb systems are plotted against lattice spacing squared 
in Figs. |3], f|, and |5] respectively. For each system the splittings obtained with ii ,L and «o,p, 
and at 0{y A ) and 0(t> 6 ), are plotted together in Figs. |3|-f)[ We collect all results with m Z/ 
in Fig. [], and all results with w 0j p m Fig- 0- 

There are a number of very clear features in the data. To begin with, we note that 
the results with mo,l for the three quarkonium systems show much smaller scaling violations 
than the results with iio,p (compare Figs. ^ and 0). The smallest scaling violations are in the 
results with u 0i l at 0(v 6 ), which show remarkably little change as the lattice spacing is varied 
here by a factor of about 2.5. The bb data show the largest scaling violations, as expected 
from the fact that this system should have the smallest size of the three. This scaling 
analysis provides evidence that mo,l tadpole renormalization yields a more continuum- like 
action than does wo,p- The improved scaling may also demonstrate indirectly that 0(v 6 ) 
corrections improve the matching of NRQCD to true QCD, as expected. 

Perhaps the most striking feature of the data is the sharp drop in the be splitting at 
smaller lattice spacings, when w 0i p is used at 0(v 6 ) (see the filled circles in Fig. [|). In fact, 
most of the c-quark data with uo,p show very large changes at the smallest lattice spacings. 
The uo,l data on the other hand exhibit a much smoother behavior. Note for example that 
the 0(f 6 ) cc splitting with m 0j p lies well below the data with u 0) l, except at the smallest 
lattice spacings, where the m 0j p data show a sharp upturn (compare the filled circles and 
filled squares in Fig. [|). 

We interpret these features as possible indicators of a breakdown in the NRQCD effective 
action at smaller lattice spacings, when the bare quark mass in lattice units clMq falls 
below about one. From Tables [III] and |V| we see that the bare c-quark mass for example 
turns out to be much larger when w 0ji is used, compared to when m ,p is used on lattices 
with comparable spacings (compare aM c ° = 0.65 at ap « 0.14 fm with aM? = 1.10 at 
cil ~ 0.16 fm). Although there is clearly a pathological behavior in the c-quark simulations 
with uo,p at smaller lattice spacings, we cannot reach definitive conclusions regarding the 
breakdown of the effective action without doing uo,l simulations at much smaller aM c °. This 
requires lattices with much smaller spacings than we can realistically simulate. 

Another key feature of these results is that the relativistic corrections to the hyperfine 
splittings are smaller when the action is renormalized using «o,l, particularly for the cc and 
be systems. For example, we find that the charmonium hyperfine splitting is reduced by 
about 30-40% in going from 0(v 4 ) to 0{y & ) when using w ,l, compared to a reduction of 
about 40-60% when using wo,p. Note that the wo,p estimate of the relativistic corrections 
depends very strongly on the lattice spacing, increasing rapidly as a decreases; this may be 
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related to pathologies in the data at oM c ° < 1, discussed above. 

The size of the corrrection to the hyperfine splitting obtained with m 0i l is consistent with 
a naive estimate of v^ 5 from operator expectation values in NRQCD ||. Hence while the 
velocity expansion for charmonium is subject to large corrections, it may not be as unreliable 
as was suggested in Ref. ||, based on simulations with uo,p on fine lattices. 

We note finally that it is reasonable to attempt to extrapolate the 0(v 6 ) hyperfine 
splittings for cc and be to zero lattice spacing, from the data on coarse lattices, where there 
is reasonably good scaling behavior (and where cutoff effects in the effective theory should 
be small, since clMq > 1 in this region However, the extrapolations in the Mo,l and 

Uq^p data, both of which exhibit good scaling on coarse lattices, are clearly very different 
(compare the filled circles with the filled squares in Fig. |], an d in Fig. £|). This suggests 
that some relevant operator coefficients q in the NRQCD action (Eqs. (|8]) and (||)) receive 
significant 0(a s ) corrections in one or both of the two tadpole schemes. This underlines the 
need to go beyond tree-level tadpole improvement in order to fully clarify the differences 
between renormalization schemes. 

IV. SUMMARY 

We have presented new evidence that clearly favors tadpole renormalization using the 
mean-link in Landau gauge over the fourth root of the average plaquette. This includes a 
demonstration of much better scaling behavior of the hyperfine splittings in three quarko- 
nium systems when uq,l is used, and a smaller size for spin-dependent relativistic corrections. 
The results presented here also help to elucidate the structure of the NRQCD effective ac- 
tion. In particular, we see signs of a breakdown in the NRQCD expansion when the bare 
quark mass falls below about one in lattice units, with pathological behavior clearly visible in 
the c-quark systems with w 0i p tadpoles. With m 0j l on the other hand the bare quark masses 
turn out to be much larger than with «o,p, resulting in a much smoother behavior. We 
have in fact been unable to do realistic simulations on lattice fine enough to make aM c ° < 1 
with m 0j l tadpoles. At the same time, these results also serve to underline the need to go 
beyond tree-level matching of improved actions, in order to eliminate uncertainties due to 
uncalculated 0(a s ) renormalizations fl8| . |l9fl . 
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TABLES 



fh (gReTr U^) (|ReTr Z7 pl ) 1/4 a eg (fm) Volume 

7.5 0.836 0.879 0.16 12 3 x 16 
7.4 0.829 0.875 0.18 10 3 x 16 
7.0 0.780 0.850 0.28 6 3 x 10 
6.85 0.763 0.840 0.32 6 3 x 10 
6.7 0.750 0.830 0.36 6 3 x 10 

6.6 0.743 0.825 0.38 6 3 x 10 



TABLE I. Simulation parameters using the Landau gauge mean-link to determine the tadpole 
renormalization uo,l (second column). The lattice spacing determined from the charmonium system 
is given as a guide. 



Pp (gReTr U^) (|ReTr Z7 pl ) 1/4 a eg (fm) Volume 

7.3 0.837 0.878 0.14 12 3 x 16 
7.2 0.830 0.875 0.17 10 3 x 16 
7.0 0.812 0.864 0.21 8 3 x 10 
6.8 0.791 0.854 0.26 6 3 x 10 
6.6 0.771 0.841 0.3 1 6 3 x 10 

6.4 0.753 0.829 0.35 6 3 x 10 
6.25 0.741 0.821 0.39 6 3 x 10 



TABLE II. Simulation parameters using the average plaquette to determine the tadpole renor- 
malization uq^p (third column). 



cc be bb 





a (fm) 


aM c °[n] 


a (fm) 




a (fm) 


aM b °[n] 


7.5 


0.155(4) 


1.10 [4] 


0.138(4) 


3.20[2] 


0.128(3) 


3.20[2] 


7.4 


0.179(2) 


1.20 [4] 


0.161(2) 


3.57[2] 


0.152(2) 


3.57[2] 


7.0 


0.280(4) 


1.97[2] 


0.261(4) 


6.10[2] 


0.257(3) 


5.35[2] 


6.85 


0.319(5) 


2.25[2] 


0.299(5) 


6.50[2] 


0.296(4) 


5.90[2] 


6.7 


0.361(6) 


2.50[2] 


0.339(6) 


7.20[2] 


0.343(6) 


6.35[2] 


6.6 


0.380(7) 


2.67[2] 


0.363(7) 


7.50[2] 


0.363(7) 


6.66[2] 



TABLE III. Lattice spacings and bare quark masses for the three quarkonium systems at 0(i> 6 ), 
using Landau gauge mean- link tadpoles uo,l; the stability parameter n for each mass is given in 
square brackets. 
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cc be bb 



Pp 


a (fm) 


aM c °[n] 


a (fm) 


aM°[n] 


a (fm) 


aM°[n] 


7.3 


0.140(4) 


0.65 [8] 


0.131(4) 


2.87[2] 


0.127(4) 


2.87[2] 


7.2 


0.169(2) 


0.83[4] 


0.150(2) 


3.20[2] 


0.145(2) 


3.20[2] 


7.0 


0.210(2) 


1.10 [4] 


0.191(2) 


4.10[2] 


0.185(2) 


3.95[2] 


6.8 


0.256(3) 


1.43 [3] 


0.235(3) 


4.98[2] 


0.228(3) 


4.53[2] 


6.6 


0.313(4) 


1.80 [3] 


0.288(4) 


5.83[2] 


0.284(4) 


5.23 [2] 


6.4 


0.350(6) 


2.15[2] 


0.328(6) 


6.45 [2] 


0.328(6) 


5.60[2] 


6.25 


0.390(6) 


2.41[2] 


0.363(7) 


6.85[2] 


0.362(6) 


5.99[2] 



TABLE IV. Lattice spacings and bare quark masses for the three quarkonium systems at 0(v 6 ) 
using average plauqette tadpoles uo,p. 



t ■ It 

nun / iiicix 


'Pi 


3 5x 


^0 


3 Si - 1 S 


2/16 


0.627(9) 


0.306(2) 


0.285(1) 


0.0210(3) 


3/16 


0.621(10) 


0.305(2) 


0.284(1) 


0.0212(3) 


4/16 


0.618(12) 


0.305(2) 


0.284(1) 


0.0212(4) 


5/16 


0.615(14) 


0.305(2) 


0.284(2) 


0.0210(4) 


6/16 


0.612(17) 


0.304(2) 


0.283(2) 


0.0209(4) 


7/16 


0.606(20) 


0.304(2) 


0.283(2) 


0.0207(5) 


8/16 


0.598(24) 


0.303(2) 


0.283(2) 


0.0205(5) 


9/16 


0.586(30) 


0.303(2) 


0.283(2) 


0.0205(6) 


10/16 


0.591(39) 


0.304(2) 


0.283(2) 


0.0204(7) 




TABLE V. Examples of fits to the 0(v 6 ) be spectra at L = 7.5 (a « 


.14 fm). 


^min/^max 


'Pi 




'So 


3 ^i " ^0 


2/16 


0.879(9) 


0.578(2) 


0.581(2) 


-0.0029(3) 


3/16 


0.874(10) 


0.577(2) 


0.579(2) 


-0.0028(3) 


4/16 


0.872(12) 


0.577(2) 


0.579(2) 


-0.0027(3) 


5/16 


0.870(14) 


0.577(2) 


0.580(2) 


-0.0026(4) 


6/16 


0.871(17) 


0.577(2) 


0.580(2) 


-0.0026(4) 


7/16 


0.877(21) 


0.579(2) 


0.581(2) 


-0.0026(5) 


8/16 


0.885(25) 


0.579(2) 


0.582(2) 


-0.0027(5) 



TABLE VI. Examples of fits to the 0(v 6 ) be spectra at (3 P = 7.3 (a w .13 fm). 
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Pl 


'Pi 




3 s 1 


1 S'o 




3 c 1 Q 
Di — Do 


7.5 

7.4 

7.0 

6.85 

6.7 

6.6 


0.62(1) 
0.70(1) 
0.95(1) 
1.01(1) 
1.08(2) 
1.12(2) 




305(2) 
0.330(1) 
0.346(1) 
0.329(1) 
0.308(1) 
0.292(1) 


284(1) 
0.306(1) 
0.310(1) 
0.287(1) 
0.261(1) 
0.245(1) 




0212(3) 
0.0235(2) 
0.0356(3) 
0.0425(3) 
0.0460(4) 
0.0471(3) 


TABLE VII 


. Final fit results for the dimensionless ener^ 


nes for the 0(v 6 


) be spectra using uq ; l- 


Pp 


'Pi 




3 S 1 


'So 




3 S, - 'So 


7.3 
7.2 
7.0 
6.8 
6.6 
6.4 
6.25 


0.87(1) 
0.98(1) 
1.12(1) 
1.25(1) 
1.36(1) 
1.44(2) 
1.48(2) 




0.577(2) 
0.637(1) 
0.681(1) 
0.706(1) 
0.699(1) 
0.681(1) 
0.648(1) 


0.580(2) 
629(1) 
0.666(1) 
0.687(1) 
0.676(1) 
0.654(1) 
0.621(1) 




-0.0027(3) 
0.0074(1) 
0.0138(1) 
0.0189(2) 
0.0233(2) 
0.0263(2) 
0.0277(2) 


TABLE VIII. Final fit results for the dimensionless energies for the 0(v 6 


) be spectra using uq,p- 


Pl 


cc 


0(v 6 ) 
be 


bb 


cc 


0(t> 4 ) 
be 


bb 


7.5 

7.4 

7.0 

6.85 

6.7 

6.6 


56.0(24) 

53.8(9) 

52.0(11) 

52.4(12) 

50.1(12) 

48.9(11) 


30.2(14) 

28.7(5) 

26.8(6) 

28.0(7) 

26.8(7) 

25.6(6) 


28.9(12) 

26.8(4) 

23.2(5) 

22.0(5) 

20.3(5) 

19.2(5) 


84.5(15) 
76.9(17) 

69.8(17) 


47.0(9) 
38.3(9) 

34.0(9) 


31.9(5) 
25.0(5) 

21.5(5) 



TABLE IX. Final fit results for hyperfine splittings in MeV for the three quarkonium systems, 
with uo,l tadpole renormalization. 
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pp 


cc 


U(v ) 
be 


00 


cc 


U(y ) 
be 


00 


Y.o 


76.4(36) 


-4.2(3) 


22.7(10) 








7.2 


39.8(7) 


9.7(2) 


21.5(4) 


100.4(18) 


45.2(8) 


26.3(4) 


7.0 


34.9(6) 


14.2(3) 


19.2(4) 








6.8 


33.0(7) 


15.8(3) 


17.1(3) 


72.8(15) 


32.1(7) 


19.4(4) 


6.6 


30.7(6) 


16.0(4) 


14.6(3) 








6.4 


29.3(7) 


15.8(4) 


13.3(3) 








6.25 


27.7(6) 


15.1(4) 


12.2(3) 


46.8(11) 


22.3(6) 


13.8(3) 



TABLE X. Final fit results for hyperfine splittings in MeV for the three quarkonium systems, 
with uq,p tadpole renormalization. 
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FIGURES 

FIG. 1. Effective mass plot for 0(v 6 ) be spectra at (3 L = 7.5 (a ss 0.14 fm): (a) 1 P 1 state (□) 
and 1 So state (o); (b) hyperfine splitting. 

FIG. 2. Effective mass plot for 0(v 6 ) be spectra at P = 7.3 (a « 0.13 fm): (a) x Pi state (□) 
and 1 S'o state (o); (b) hyperfine splitting. 



FIG. 3. Hyperfine splittings for the cc system versus lattice spacing squared. 
FIG. 4. Hyperfine splittings for the be system versus lattice spacing squared. 
FIG. 5. Hyperfine splittings for the bb system versus lattice spacing squared. 

FIG. 6. Hyperfine splittings with uq^l versus lattice spacing squared. 

FIG. 7. Hyperfine splittings with uq,p versus lattice spacing squared. 
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Figure 1b 
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Figure 2a 
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Figure 2b 
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